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Abstract 

We extend the wellposcdncss results for second order backward stochastic differential 
equations introduced by Soner, Touzi and Zhang |24j to the case of a bounded terminal 
condition and a generator with quadratic growth in the z variable. More precisely, 
we obtain uniqueness through a representation of the solution inspired by stochastic 
control theory, and we obtain two existence results using two different methods. In 
particular, we obtain the existence of the simplest purely quadratic 2BSDEs through 
the classical exponential change, which allows us to introduce a quasi-sure version of 
the entropic risk measure. As application, we also study robust risk-sensitive control 
problems. Finally, we prove a Feynman-Kac formula and a probabilistic representation 
for fully nonlinear PDEs in this setting. 



Key words: Second order backward stochastic differential equation, quadratic growth, 
r.c.p.d., Feyman-Kac, fully nonlinear PDEs, quasi-sure. 

AMS 2000 subject classifications: 60H10, 60H30 



'Research supported by the Chair Financial Risks of the Risk Foundation sponsored by Societe Generale, 
the Chair Derivatives of the Future sponsored by the Federation Bancaire Francaise, and the Chair Finance 
and Sustainable Development sponsored by EDF and Calyon. 

^CMAP, Ecole Poly technique, Paris, dylan.possamai@polytechnique.edu. 

*CMAP, Ecole Polytechnique, Paris, chao.zhou@polytechnique.edu. 



1 



1 Introduction 

Backward stochastic differential equations (BSDEs for short) appeared for the first time in 
Bismut [3] in the linear case. However, they only became a popular field of research after 
the seminal paper of Pardoux and Peng [20], mainly because of the very large scope of their 
domain of applications, ranging from stochastic control to mathematical finance. 

On a filtered probability space (Q,,J-,{J-t} 0<t<T ,¥) generated by an K^-valued Brownian 
motion B, solving a BSDE amounts to finding a pair of progressively measurable processes 
(Y, Z) such that 



where / (called the generator) is a progressively measurable function and £ is an Tt- 
measurable random variable. 

Pardoux and Peng [20] proved existence and uniqueness of the above BSDE provided that 
the function / is uniformly Lipschitz in y and z and that £ and / s (0, 0) are square integrable. 
Then, in [21] , they proved that if the randomness in / and £ is induced by the current value 
of a state process defined by a forward stochastic differential equation, then the solution to 
the BSDE could be linked to the solution of a semilinear PDE by means of a generalized 
Feynman-Kac formula. This link then opened the way to probabilistic numerical methods 
for solving semilinear PDEs, which particularly well suited for highly dimensional problems 
(see Bouchard and Touzi [3] among many other). 

Nonetheless, the class of fully nonlinear PDEs remained inaccessible when considering only 
BSDEs, until a recent work of Cheredito, Soner, Touzi and Victoir [9]. They introduced a 
notion of second order BSDEs (2BSDEs), which were then proved to be naturally linked to 
fully nonlinear PDEs. However, only a uniqueness result in the Markovian case was proved, 
and no existence result (apart from trivial ones) were available. Following this work, Soner, 
Touzi and Zhang [21] gave a new definition of 2BSDEs, and thus provided a complete 
theory of existence and uniqueness under uniform Lipschitz conditions similar to those of 
Pardoux and Peng. Their key idea was to reinforce the condition that the 2BSDE must 
hold P — a.s. for every probability measure P in a non-dominated class of mutally singular 
measures (see Section [2] for precise definitions). Let us describe briefly the intuition behind 
their definition. 

Suppose that we want to study the following fully non-linear PDE 



If the function 7 1— > h(t,x,r,p,j) is assumed to be convex, then it is equal to its double 
Fenchel-Legendre transform, and if we denote its Fenchel-Legendre transform by /, we have 




du 
~dt 



h (t, x, u(t, x),Du(t, x), D 2 u(t, 2)) = 



u(T,x) = g(x). 



(1.1) 




(1.2) 
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Then, from (jl.2p . we expect, at least formally, that the solution u of (jl.ip is going to verify 



u(t,x) = sup u a (t,x) 

a>0 



where u a 



is defined as the solution of the following semi-linear PDE 




u a (T,x)=g(x). (1.3) 



Since u a is linked to a classical BSDE, the 2BSDE associated to u should correspond (in 
some sense) to the supremum of the family of BSDEs indexed by a. Furthermore, changing 
the process a can be achieved by changing the probability measure under which the BSDE 
is written. However, this amounts to changing the quadratic variation of the martingale 
driving the BSDE, and therefore leads to a family of mutually singular probability mea- 
sures. In these respects, the 2BSDE theory shares deep links with the theory of quasi-sure 
stochastic analysis of Denis and Martini p.0] and the theory of G-expectation of Peng [22]. 

Following the breakthrough of [23], Possamai [23] extended the existence and uniqueness 
result for 2BSDEs to the case of a generator having linear growth and satisfying a mono- 
tonicity condition. Motivated by a robust utility maximization problem under volatility 
uncertainty (see the accompanying paper [19]), our aim here is to go beyond the results of 
[23 j to prove an existence and uniqueness result for 2BSDEs whose generator has quadratic 
growth in z. 

The question of existence and uniqueness of solutions to these quadratic equations in the 
classical case was first examined by Kobylanski [T7], who proved existence and uniqueness 
of a solution by means of approximation techniques borrowed from the PDE litterature, 
when the generator is continuous and has quadratic growth in z and the terminal condition 
£ is bounded. Then, Tevzadze [28] has given a direct proof for the existence and uniqueness 
of a bounded solution in the Lipschitz-quadratic case, proving the convergence of the usual 
Picard iteration. Following those works, Briand and Hu [6] have extended the existence 
result to unbounded terminal condition with exponential moments and proved uniqueness 
for a convex coefficient [7] . Finally, Barrieu and El Karoui [2] recently adopted a completely 
different approach, embracing a forward point of view to prove existence under conditions 
similar to those of Briand and Hu. 

In this paper, we propose two very different methods to prove the wellposedness in the 
2BSDE case. First, we recall some notations in Section [2] and prove a uniqueness result 
in Section [3] by means of a priori estimates and a representation of the solution inspired 
by the stochastic control theory. Then, Section HJ is devoted to the study of approxima- 
tion techniques for the problem of existence of a solution. We advocate that since we are 
working under a family of non-dominated probability measures, the monotone or domi- 
nated convergence theorem may fail. This is a major problem, and we spend some time 
explaining why, in general, the classical methods using exponential changes fail for 2BSDEs. 
Nonetheless, using very recent results of Briand and Elie [8], we are able to show a first 
existence result using an approximation method. Then in Section [5] we use a completely 
different method introduced by Soner, Touzi and Zhang [24] to construct the solution to 
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the quadratic 2BSDE path by path. Next, we use these results in Section [6] to study an 
application of 2BSDEs with quadratic growth to robust risk-sensitive control problems. Fi- 
nally, in Section [71 we extend the results of Soner, Touzi and Zhang |24j on the connections 
between fully non-linear PDEs and 2BSDEs to the quadratic case. 

2 Preliminaries 

Let !J:= {w£ C([0,T],R d ) : lu = 0} be the canonical space equipped with the uniform 
norm ||w|| := sup 0<t<T \cot\, B the canonical process, Po the Wiener measure, F := 
{-^i}o<t<T * ne filtration generated by B, and F + := {-T 7 ^} 0<t<T the right limit of F. 

2.1 A first set of probability measures 

Our aim here is to give a correct mathematical basis to the intuitions we provided in the 
Introduction. We first recall that by the results of Karandikar |15j . we can give pathwise 
definitions of the quadratic variation (B) t and its density at- 

Let Vw denote the set of all local martingale measures P (i.e. the probability measures P 
under which B is a local martingale) such that 

(B) t is absolutely continuous in t and a takes values in S^ -0 , P — a.s. (2.1) 

where S^ -0 denotes the space of all d x d real valued positive definite matrices. ). 

As in [23], we concentrate on the subclass Vs C Vw consisting of all probability measures 

P° := P o (X")- 1 where X t ° := f a] /2 dB s , t € [0,T], P - a.s. (2.2) 

Jo 

for some F-progressively measurable process a satisfying J^\a s \ds < +oo. We recall 
from [25] that every P € Vs satisfies the Blumenthal zero-one law and the martingale 
representation property. 

Notice that the set Vs is bigger than the set Vs introduced in [23], which is defined by 

V s := {P Q G V s , a < a < a, P - a.s.} , (2.3) 
for fixed matrices a and a in S^ -0 . 

2.2 The Generator and the final set Vh 

Before defining the spaces under which we will be working or defining the 2BSDE itself, we 
first need to restrict one more time our set of probability measures, using explicitely the 
generator of the 2BSDE. 

Following the PDE intuition recalled in the Introduction, let us first consider a map 
H t {uj,y,z,i) ■ [0,71 x n x R x R d x D H -> R, where D H C R dxd is a given subset 
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containing 0. As expected, we define its Fenchel-Legendre conjugate w.r.t.7 by 

F t (u,y,z,a) := sup j -Tr(o7) - H t (u, y, z, 7) 1 for a G 
ieD H J 

F t (y,z) := F t (y,z,a t ) and F? := F t (0,0). 

We denote by D Ft ^ y z ^ the domain of F in a for a fixed (i, u>, y, z), and as in [23] we restrict 
the probability measures in Vu C "Ps 

Definition 2.1. Pjj consists of all ¥ 6 suc/t i/iat 

a P < a < ap, dt x dF — a.s. for some a P , ap G S^°, and € D Ft r y>z y 

Remark 2.1. XTie restriction to the setVu obeys two imperatives. First, since F is destined 
to be the generator of our 2BSDE, we obviously need to restrict ourselves to probability 
measures such that at € D Ft r y>z \. Moreover, we also restrict the measures considered to the 
ones such that the density of the quadratic variation of B is bounded to ensure that B is 
actually a true martingale under each of those probability measures. This will be important 
to obtain a priori estimates. 

2.3 Spaces of interest 

We now provide the definition of the spaces which will be used throughout the paper. 

2.3.1 Quasi-sure spaces 

denotes the space of all Fr-measurable scalar r.v. £ with 

U\k°s := SU P ll£IL~(P) < +°°- 



M P H denotes the space of all F + -progressively measurable M^-valued processes Z with 

< +00. 



\Z\f nP := supE* 



VV2 7 |2 ^ 1 

\a t Zt\ dt 







denotes the space of all F + -progressively measurable R-valued processes Y with 
Vh — cadlag paths, and H^H^oo := sup H^H^oo < +00. 

H 0<t<T H 

Finally, we denote by UCfe(Q) the collection of all bounded and uniformly continuous maps 
£ : £1 — > K with respect to the ||-|| -norm, and we let 

Clfi := the closure of UCfe(ri) under the norm ||-|| L ^. 

Remark 2.2. We emphasize that all the above norms and spaces are the natural general- 
izations of the standard ones when working under a single probability measure. 
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2.3.2 The space BMO(?h) and important properties 



It is a well known fact that the Z component of the solution of a quadratic BSDE with 
a bounded terminal condition belongs to the so-called BMO space. Since this link will 
be extended and used intensively throughout the paper, we will recall some results and 
definitions for the BMO space, and then extend them to our quasi-sure framework. We 
first recall (with a slight abuse of notation) the definition of the BMO space for a given 
probability measure P. 

Definition 2.2. BMO(P) denotes the space of all F + -progressively measurable W 1, -valued 
processes Z with 

cT 



\z\ 



BMO( 



sup 



E_ 



J \a] /2 Z t \ 2 dt 



< +00, 



where Tq is the set of Tt stopping times taking their values in [0, T]. 

We also recall the so called energy inequalities (see |16] and the references therein). Let 
Z € BMO(P) and p > 1. Then we have 



T 



iV 2 z f 



ds 



< 2p\ [A\\Z 



BMO 



(2.4) 



The extension to a quasi-sure framework is then naturally given by the following space. 
BMO('Ph) denotes the space of all F + -progressively measurable Revalued processes Z with 



Ibmo(Ph) 



SU P II-^IIbmo(p) < +°°- 



The main interest of the BMO spaces is that if a process Z belongs to it, then the stochastic 
integral J Q Z s dB s is a uniformly integrable martingale, which in turn allows us to use it for 
changing the probability measure considered via Girsanov's Theorem. The two following 
results give more detailed results in terms of L r integrability of the corresponding Doleans- 
Dade exponentials. 

Lemma 2.1. Let Z G BMO(T'h). Then there exists r > 1, such that 



sup E 

¥eV H 



£ 



Z H dB< 



< +00. 



Proof. By Theorem 3.1 in [16] , we know that if ||^||bmo(p) — ^K r ) f° r some one-to-one 
function $ from (1, +00) to R* + , then £ (f Q Z s dB s ) is in Z/(P). Here, since Z G BMO(T'h), 
the same r can be used for all the probability measures. □ 

Lemma 2.2. Let Z € BMO(P H )- Then there exists r > 1, such that for all t € [0, T] 

£ (j*Z s dB s 



sup Ef 



£ 



T 



ZAB 




< +00. 
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Proof. This is a direct application of Theorem 2.4 in |16| for all P G Vh- ^ 

We emphasize that the two previous Lemmas are absolutely crucial to our proof of unique- 
ness and existence. Besides, they also play a major role in our accompanying paper [19j . 

2.4 The definition of the 2BSDE 

Everything is now ready to define the solution of a 2BSDE. We shall consider the following 
2BSDE 

Yt = Z + J F S (Y S , Z s )ds - J Z s dB s + K T -K t , 0<t<T, V H - qs. (2.5) 

Definition 2.3. We say (Y, Z) G H<g x M H is a solution to 2BSDE (J23]) if 

• Y T = i,Vu - q.s. 

• For all P G Vh, the process K ¥ defined below has non- decreasing paths P — a.s. 

Kf:=Y -Y t - [ F s {Y s ,Z s )ds + [ Z s dB s , < t < T, P — a.s. (2.6) 
Jo Jo 

• The family jf£T p ,P G Vh} satisfies the minimum condition 

Kf = essinf p Ef' \k£] , < t < T, P - a.s., VP G V H , (2.7) 

P H (i + ,IP) := {f' € Ph, P' =P on 7^+} . 

Moreover if the family {K ¥ ,F G "Ph} can 6e aggregated into a universal process K , we call 
(Y, Z, K) a solution of 2BSDE ([23]) . 

Remark 2.3. Let us comment on this definition. As already explained, the PDE intuition 
leads us to think that the solution of a 2BSDE should be a supremum of solution of standard 
BSDEs. Therefore for each P, the role of the non- decreasing process is in some sense 
to "push " the process Y to remain above the solution of the BSDE with terminal condition 
^ and generator F under P. In this regard, 2BSDEs share some similarities with reflected 
BSDEs. 



Pursuing this analogy, the minimum condition (|2.7p tells us that the processes K act in 
a "minimal" way (exactly as implied by the Skorohod condition for reflected BSDEs), and 
we will see in the next Section that it implies uniqueness of the solution. Besides, if the set 
Vh wa s reduced to a singleton {P}, then (|2.7|) would imply that K p is a martingale and a 
non- decreasing process and is therefore null. Thus we recover the standard BSDE theory. 



Finally, we would like to emphasize that in the language of G-expectation of Peng 
\2. ?] ) is equivalent, at least if the family can be aggregated into a process K , to saying that 
—K is a G- martingale. This link has already observed in \26$ where the authors proved 
the G-martingale representation property, which formally corresponds to a 2BSDE with a 
generator equal to 0. 
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2.5 Assumptions 

We finish this Section by giving the properties which will be assumed to be satisfied by the 
generator F of the 2BSDE. 

Assumption 2.1. (i) The domain Dp t ( VtZ \ = Dp t is independent of (u,y, z). 
(ii) For fixed (y, z, 7), F is ¥ -progressively measurable. 
(hi) F is uniformly continuous in co for the \ \ ■ ||oo norm. 

(iv) F is continuous in z and has the following growth property. There exists (a, f3, 7) € 
M + x M + x R* such that 
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\F t (uj,y,z,a)\ <a + /3\y\ + - 



a l ' 2 z 



a l / 2 z 



for all (t, y, z, oj, a). 

(v) F is C 1 in y and C 2 in z, and there are constants r and 9 such that for all (t, co, y, z, a), 
\D y F t {uj,y,z,a)\ < r, \D z F t (u,y, z,a)\ <r + 6 
\D 2 zz F t (u,y,z,a)\<6. 

Remark 2.4. Let us comment on the above assumptions. Assumptions \2.1\ (i) and (hi) 
are taken from \21$ and are needed to deal with the technicalities induced by the quasi-sure 
framework. Assumptions \2. 1\ (ii) and (iv) are quite standard in the classical BSDE littera- 
ture. Finally, Assumption \2. 1\ (v) was introduced by Tevzadze in 128)1 for quadratic BSDEs. 
It allowed him to prove existence of quadratic BSDEs through fixed point arguments. This 
is this consequence which will be used for technical reasons in Section^ 

However, it was also proved in \28§ . that if both the terminal condition and F° were small 
enough, then Assumption \2. 1\ (v) can be replaced by a weaker one. We will therefore some- 
times consider 

Assumption 2.2. Let points (i) to (iv) of Assumption [2U\ hold and 

(v) We have the following "local Lipschitz" assumption in z, 3fi > and a progressively 
measurable process (j) £ BMO('Ph) such that for all (t,y, z, z',ui, a), 



F t (u,y,z,a) - F t (u,y,z',a) - (p t -a 1/2 (z - z') 



< va 1 ' 2 



z — z 



a^z 



+ 



(vi) We have the following uniform Lipschitz-type property in y 

\F t (u,y,z,a) - F t {u , y' , z , a)\ < C\y-y'\ , for all (y,y',z,t,u,a). 

Furthermore, we observe that our subsequent proof of uniqueness of a solution of our 
quadratic 2BSDE only uses Assumption (j2.2f) . 
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Remark 2.5. By Assumption \2. H v) . we have that F® is actually bounded, so the strong 
integrability condition 




< +00, 



with a constant k € (1,2] introduced in \2J$ is not needed here. Moreover, this implies that 
a always belongs to Dp t and thus that Vh *s n °t empty. 



3 A priori estimates and uniqueness of the solution 



Before proving some a priori estimates for the solution of the 2BSDE (|2.5p . we will first 
prove rigorously the intuition given in the Introduction saying that the solution of the 
2BSDE should be, in some sense, a supremum of solution of standard BSDEs. Hence, 
for any P 6 Vh, F-stopping time r, and J 7 ,— measurable random variable £ G L°°(P), we 
define (y ¥ , z ¥ ) := (y p (r, £), z ¥ (t, £)) as the unique solution of the following standard BSDE 
(existence and uniqueness have been proved under our assumptions by Tevzadze in |28j ) 

Vt = £ + f T F-M, z^)ds - f T z^dB s , < t < r, P - a.s. (3.1) 
Jt Jt 

First, we introduce the following simple generalization of the comparison Theorem proved 
in |28j (see Theorem 2). 

Proposition 3.1. Let Assumptions \2.I\ hold true. Let £1 and £2 £ L°°(P) for some prob- 
ability measure P, and V 1 , i = 1,2 be two adapted, cddldg non- decreasing processes null at 
0. Let (Y^Z*) E B°°(P) x M 2 (P), i = 1,2 be the solutions of the following BSDEs 

Yl = C+ [ Fstf, Z\)ds - [ Z\dB s + Vi- Vf, P - a.s., i = 1, 2, 
Jt Jt 

respectively. If £1 > £2, P — a.s. and V 1 — V 2 is non- decreasing, then it holds P — a.s. that 
for all t € [0, T] 

Yt 1 > Yl 

Proof. First of all, we need to justify the existence of the solutions to those BSDEs. 
Actually, this is a simple consequence of the existence results of Tevzadze |28| and for 
instance Proposition 3.1 in |18j . Then, the above comparison is a mere generalization of 
Theorem 2 in [28]. □ 
We then have similarly as in Theorem 4.4 of |24j the following results which justifies the 
PDE intuition given in the Introduction. 

Theorem 3.1. Let Assumptions \2.2\ hold. Assume £ € and that (Y, Z) G x M H is 
a solution to 2BSDE fl23|). Then, for any P G V H and0<h<t 2 < T, 

Y tl = ess sup p y({t 2 ,Y t2 ), P-a.s. (3.2) 

P'eP H (ii,P) 

Consequentlty, the 2BSDE (|2.5p has at most one solution in x M H . 
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Before proceeding with the proof, we will need the following Lemma which shows that in 
our 2BSDE framework, we still have a deep link between quadratic growth in z of the 
generator and the BMO spaces. 



Lemma 3.1. Let Assumption^. 1 hold. Assume £ G and that (Y, Z) G B'g x M 2 H 
solution to 2BSDE ([23]). Then Z G BMO(P H )- 



is a 



Proof. Denote Tq the collection of stopping times taking values in [0, T]. By Ito's formula 
under P applied to e~ vYt , which is a cadlag process, for some v > 0, we have for every r G T^f 



2 r T 



-vY t 



-1/2 ry 

a J Z t 



dt = e~< - <r vY - - v 



T j-T 



e~ vY t~dK r t -v \ e-» Yt Ft(Yt,Z t )dt 



+ v 



e~ vY t-Z t dBt 



r<S<T 



+ vAY s e~ uY s 



Since Y G B^ 1 , K p is non-decreasing and since the contribution of the jumps is negative 
because of the convexity of the function x — > e~ vx , we obtain with Assumption I2.2( iv) 



r r T 


-1/2 „ 

a t ' Z t 


2 


[I 


dt 









<e"""°» il + vT[a + P\\Y 



+ T E ' 



r r T 


-1/2 „ 

<h z t 


2 


[I ^ 


dt 









By choosing v = 27, we then have 

e! 



[/> w 


-1/2 „ 


2 
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D f (l + 2 7 T (a + /3||r|| D „)). 



Finally, we obtain 

e! 



T -1/2 5 



dt 



< 1 e 4Tl|y|| Df ( 1 + 27r ( a + / 3||y|| Ds) 



which provides the result by arbitrariness of P and r. 



□ 



Proof of Theorem 13.11 The proof follows the lines of the proof of Theorem 4.4 in 
but we have to deal with some specific difficulties due to our quadratic growth assumption. 
First if (|3.2p holds, this implies that 



Y t = esssup P yf (T, £), t G [0,T], P - a.s. for all P G Vh, 

F'£V H (t+,P) 

and thus is unique. Then, since we have that d{Y,B) t = Ztd(B) t , Vh — Q-s., Z is also 
unique. We now prove (13. 2|) in three steps. Roughly speaking, we will obtain one inequality 
using the comparison theorem, and the second one by using the minimal condition ([27 
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(i) Fix < h < t 2 < T and P G V H - For any P G V H (tf,F), we have 

Y t = Y t2 + F S (Y S , Z s )ds - Z s dB s + k( - Kf' , h<t<t 2 , P' 



a.s. 



and that K is nondecreasing, P — a.s. Then, we can apply the comparison Theorem 13.11 
under P' to obtain Y h > y£ (t 2 ,Y t2 ), P' - a.s. Since P' = P on 7"+, we get Y h > y£ (t 2 ,Y t2 ), 
P — a.s. and thus 

Y tl > esssup p y((t 2 ,Y t2 ), P — a.s. 

(ii) We now prove the reverse inequality. Fix P £ Vu- Let us assume for the time being 
that 



Cff:= esssup p E* 



K-<f 



< +oo, P — a.s., for all p > 1. 



For every P' G V H (t + ,lP), denote 

:= y - /' (i 2 , y i2 ) and 5Z := Z - / {t 2 ,Y t2 ). 

By the Lipschitz Assumption I2.2( vi) and the local Lipschitz Assumption 12. 2( v). there exist 
a bounded process A and a process r\ with 



\Vt\ < V 



-1/2 ry 
OJ Z t 



+ 



-1/2 p' 



a.s. 



such that 
SY t -- 

Define for t\ < t < t 2 



2 '\ S 5Y S + (r/ s + ^ s )a l J 2 5Z s ^ ds- J* 5Z s dB s + Kf 2 ' - Kf', t < t 2 , P' - a.s. 



M t := exp 



\ q ds 



a.s. 



Now, since (j) G BMO('Ph), by Lemma 12. 11 we know that the exponential martingale 



E ( '(<!>, + Vs)K 1/2 dB l 



is a P -uniformly integrable martingale. Therefore we can define a probability measure 
which is equivalent to P , by its Radon-Nykodym derivative 



% = £[j\<Ps + Vs))a-^dB ! 



Then, by Ito's formula, we obtain, as in |24j . that 

ft 2 



5Y tl =K 



M t dKf 



< E 



ti 



sup (M t )(Kf 2 -Kl 
ti<t<t 2 



since is non-decreasing. 
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Then, since A is bounded, we have that M is also bounded and thus for every p > 1 



E, 



sup {M t y 

t!<t<t 2 



a.s. 



(3.3) 



This is now that the we will need more work than in the Lipschitz case, and that the 
BMO properties of the solution will play a crucial role. Since (v + (f>)a s ' is in BMO(^ H ), 
we know by Lemma 12.11 that there exists r > 1, independent of the probability measure 
considered, such that 



sup E ff 



£ 



T 



+ ria)a, 



-V 2 dB. 



< +oo. 



Then it follows from the Holder inequality and Bayes Theorem that 



SY tl < 



E : 



E' 



£ 



q(4> s + r] s )a s 1/2 dB< 



E 



sup M t [Kl -Kl 

t!<t<t 2 



<C(C r t ; 







) 4 '(< 


'<-<_ 



By the minimum condition (|2.7p and since P' G Vh {t + , P) is arbitrary, this ends the proof. 

(iii) It remains to show that the estimate for Cf^ p holds for p > 1. Once again, this will be 
possible because of the BMO property satisfied by the solution of the 2BSDE. By definition 
of the family {K ¥ ,F G Vh}, we have 



E 



K r - K l 

JX t2 A tl 



<C[l + \\Y\\l^ + U\\^+Kl 



"ti 



H H 
t 2 \ P 

Z t dB t 



* 2 1/2 ^ 



dt 



Thus by the energy inequalities (|2.4|) and by Burkholder-Davis-Gundy inequality, we obtain 
that 



E 



K ¥ - K ¥ 



<C (l + \\Y\\P +||£||P +\\Z 



,2p 

BMOi 



-1- II 7\\ p 



BMO 



Therefore, we have proved that 



sup E 

v'ev H {t+,¥) 



K - K 



< +oo. 



Then we can proceed exactly as in the proof of Theorem 4.4 in 



□ 



Remark 3.1. It is interesting to notice that in contrast with standard quadratic BSDEs, for 
which the only property of BMO martingales used to obtain uniqueness is the fact that their 
Doleans-Dade exponential is a uniformly integrable martingale, we need a lot more in the 
2BSDE framework. Indeed, we use extensively the energy inequalities and the existence of 
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moments for the Doleans-Dade exponential ( which is a consequence of the so called reverse 
Holder inequalities, which is a more general version of Lemma \2.1\) . Furthermore, we will 
also use the so-called Muckenhoupt condition (which corresponds to Lemma \2.2l see 116\/ 
for more details) in both our proofs of existence. This seems to be directly linked to the 
presence of the non- decreasing processes and raises the question about whether it could 
be possible to generalize the recent approach of Barrieu and El Karoui J2j, to second- order 
BSDEs. Indeed, since they no longer assume a bounded terminal condition, the Z part of 
the solution is no-longer BMO. We leave this interesting but difficult question to future 
research. 

We conclude this section by showing some a priori estimates which will be useful in the 
sequel. Notice that these estimates also imply uniqueness, but they use intensively the 
representation formula 



Theorem 3.2. Let Assumption \2.2\ hold. 

(i) Assume that £ € and that (Y, Z) G x M 2 H is a solution to 2BSDE ([23]) . Then, 
there exists a constant C such that 



Idling + II-^IIbmo(Ph) - ^ [} + H£Hl 

(K%-k*Y] <c(i + u\\i 



Vp > i, sup e; 

¥£T h , reTj 

(ii) Assume that € L 1 ^ and that (Y\Z l ) € x ~M? H is a corresponding solution to 
2BSDE (H3D, i = 1,2. Denote 5£ := £ x - £ 2 , 5Y := Y 1 - Y 2 , 5Z := Z 1 - Z 2 and 
SK V := K^' 1 — K^' 2 . Then, there exists a constant C such that 

\\sY\\^<c\m^ 

II^IIbmo(p h ) - c ll^lkg 3 (i + II^ILf + II^ILg 3 



\fp > 1, sup E 

P£V H 



sup 

0<t<T 



6K* 



<c\\t\\lg + + 



Proof, (i) By Theorem 13.11 we know that for all P € Vh an d f° r a h t G [0, T] we have 

v p' 
it = ess sup y t . 

Then by Lemma 1 in [6], we know that for all P £ Vh 

1 , / P \ , , _ / e? T - 1 



Vt 



< — log (w^ i where ip(x) := exp (^ya- — - — - + ^fe^ T x\ . 



Thus, we obtain \yf\ < a [eP — l) //3 + IICIIl^ > an< l by the representation recalled 
above, the estimate of ||V|| DOO is obvious. By the proof of Lemma 13. 11 we have now 
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Finally, we have for all r € Tq ', for all P £ Vh an d for all p > 1, by definition 

(KV-Kly=\Y T -fi- J F t {Y y ,Zt)dt + J ZtdBA . 
Therefore, by our growth Assumption I2.2l fiv) 



E 



<c i + tt- + ni£-+E? 



<c 1 + neiifcc + i|z 



,2p 

Ibmo(-Ph) 



+ IIZI 



BMO(-Ph) 



<c i + neii 



where we used again the energy inequalities and the BDG inequality. This provides the 
estimate for K p by arbitrariness of r and P. 

(ii) With the same notations and calculations as in step (ii) of the proof of Theorem 13. 1| it 
is easy to see that for all P G Vh an d for all t € [0,T], we have 



6j$ = E?[M T 5Z]<C\\5Z\ 



since M is bounded and we have (13. 3p . By Theorem l3.lt the estimate for 5Y follows. Now 
apply Ito's formula under a fixed P € Vh to \5Y\ 2 between r € T^f and T 



E 



\5Y r \ 2 + 



T 1/2 2 



dt 



<e! 



l' T 5Y t -d{5K v t 



Then, we have by Assumption I2.2f iv) and the estimates proved in (i) above 



Ez 



aJ dZ t 



dt 



< C \\5Y || D g 

+ ll^fe+2 

< C 



i + EII^IU + ll^ 



i=l 



i + II^IU + He 2 



BMO(-Ph) 

+ 



K ¥ /-K^ 2 



hi 00 / ' 



which implies the required estimate for 5Z. Finally, by definition, we have for all P € Vh 
and for all t £ [0, T] 



SK[ = 5Y - 5Y t - J F S (Y S L , Z l s ) - F S (Y S \ Z z s )ds + / 6Z s dB s . 
By Assumptions [2T2T v) and (vi), it follows that 



sup 

0<t<T 



5K 



T 

< C I \\5Y\\ 000 + / \a]6Z s \(l + \a]Z]\ + \a] Z 2 \)ds + sup 

0<t<T 



5ZedB< 



and by Cauchy-Schwarz, BDG and energy inequalities, we see that 
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E 







p 






sup 


SKf 




<CE ¥ 


u 


0<t<T 









l + \alzl\ 2 + \a\z 2 \ 2 ds 











(J \aUz s \ 2 d^j 















[or 


al /2 5Z s 







2 \ Vl 2 ' 

ds 



<c\mti{i^\e\C+\\eti 



□ 

Remark 3.2. Let us note that the proof of (i) only requires that Assumption \2.2{ \v) holds 
true, whereas (ii) also requires Assumption \2.2( v) and (vi). 

4 2BSDEs and monotone approximations 

This Section is devoted to the study of monotone approximations in the 2BSDE framework. 
We start with the simplest quadratic 2BSDEs, which allows us to introduce a quasi-sure 
version of the entropic risk measure. In that case, we obtain existence through the classical 
exponential change. Then, we show that for more general generators, this approach usually 
fails because of the absence of a general quasi-sure monotone convergence Theorem. Finally, 
we prove an existence result using another type of approximation which has the property 
to be stationary. 



4.1 Entropy and purely quadratic 2BSDEs 

Given ^ E we first consider the purely quadratic 2BSDE defined as follows 

Then we use the classical exponential change of variables and define 

f-t 



J Z s dB s + — Kf, < t < T, Vh- q*. (4.1) 



Y t := e 



Ft 

, Z t := jYtZt, Zf := 7 / Y a dK* - £ tf Y ' - e* Y .- - "fAY 6 

0<s<t 



At least formally, we see that (Y, Z, K ) verifies the following equation 

rT m 



Y t = e 



-75 



Z s dB s + K T — K t , < t < T, P - a.s. VP G V H 



(4.2) 



which is in fact a 2BSDE with generator equal to (and thus Lipschitz), provided that 



the familly ( K I satisfies the minimum condition (12.71) . Thus the purely quadratic 

2BSDE (|4.ip is linked to the 2BSDE with Lipschitz generator (|4.2h . which has a unique 
solution by Soner, Touzi and Zhang [23]. We now make this rigorous. 
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Proposition 4.1. The 2BSDE (f^TTj) has a unique solution (Y,Z) G Bg> x M. 2 H given by 



Y t = - In ( ess sup p E P ' 

T \P'e7>ff(t+,P) 



a.s., f G [0,T], /or a// P G V H - 



Proof. Uniqueness is a simple consequence of Theorem 13. 11 In the following, we prove the 
existence in 3 steps. 

Step 1: Let (Y,Z) G TD 2 H x M 2 H be the unique solutinon to the 2BSDE (@3]) and if be 
the correponding non-decreasing processes. In particular, we know that 



Yt = ess sup ssL t 



-7* 



a.s., 



which implies that Y G B 1 ^, since 

< e 



■7llfll L ~ <y t < e TllfllLf _ 



We can therefore make the following change of variables 

F t :=iln(F t ), Z t :=-^. 

7 7Y t 

Then by Ito's formula, we can verify that the pair (Y, Z) G B 1 ^ x satisfies (|4.ip with 



Iff 



o 7^ 



Moreover, notice that is non-decreasing with Kq = 0. 

Step 2: Denote now (y p ,z p ) the solutions of the standard BSDEs corresponding to the 
2BSDE (|4.ip (existence and uniqueness are ensured for example by [2Sj). Furthermore, if 
we define 



lit 



— P P 

7y* ^ , 



then we know that (y p , z p ) solve the standard BSDE under P corresponding to (|4.2p . Due 
to the monotonicity of the function x — > ln(x) and the representation for Y 



it = ess sup y t = ess sup SL t 

P'eP H (i+,P) P'GP H (i+,P) 

we have the following representation for Y 



1 



a.s., 



Yt = ess sup y t = — In ess sup E t 

P'eV H (t+,f) ^ \P'G-PH(i+,P) 



a.s. 



Step 3: Finally, it remains to check the minimum condition for the family of non-decreasing 
processes Since the purely quadratic generator satisfies the Assumption 12 . T\ we can 

derive the minimum condition from the above representation for Y exactly as in the proof 
of Theorem 14.11 in Subsection 14.31 □ 
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Thanks to the above result, we can define a quasi-sure (or robust) version of the entropic 
risk measure under volatility uncertainty 



e 7 t (£) := - In ( esssup p E P ' 



-7? 



where the parameter 7 stands for the risk tolerance. We emphasize that, as proved in [25 
(see Proposition 4.11), the solution of (|4.ip is actually F- measurable, so we also have 



e 7,t(£) := — In ( esssup P E P 
T \F'e , PH(t,P) 



-7« 



which in particular implies that 



1 



e 7 o(£) = — In sup E 
7 \^eP„ 



-75 



More generally, by the same exponential change and arguments above, we can also prove 
that there exists a unique solution to 2BSDEs with terminal condition £ 6 and the 



following type of quadratic growth generators a 1//2 z<?(i, oS) + h(t, uj) + 
h are assumed to be bounded, adapted and uniformly continuous in uj for the 



^1/2 
aJ z 



where g and 



4.2 Why the exponential transformation may fail in general ? 

Coming back to Kobylanski [T7], we know that the exponential transformation used in the 
previous subsection is an important tool in the study of quadratic BSDEs. However, unlike 
with a purely quadratic generator, in the general case the exponential change does not lead 
immediately to a Lipschitz BSDE. For the sake of clarity, let us consider the 2BSDE (12. 5j) 
and let us denote 



r, := e* Y t := e** Z t := 7^, «f := 7 / - E ^ " ^ " 7 AF * 

j „:~"7, 



0<s<< 



Then we expect that, at least formally, if (Y,Z) is a solution of (|2.5p . then (Y,Z) is a 
solution of the following 2BSDE 

2' 



Y t = V + 1 Y 



( 



V 



foe y. s z. 



7 lY s 



a s Zj s 



ds 



Z s dB s + K T - K t . (4.3) 



Let us now define for (t, y, z) £ [0, T] x 

/ 

G t {uj,y,z) := jy 

\ 



logy z 
F t uj, , — 

7 IV 



?V2 



2 7 ^ 



Then, despite the fact that the generator G is not Lipschitz, it is possible, as shown by 
Kobylanski [17], to find a sequence (G n ) n >o of Lipschitz functions which decreases to G. 
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Then, it is possible, thanks to the result of [23] to define for each n the solution (Y n , Z n ) of 
the corresponding 2BSDE. The idea is then to prove existence and uniqueness of a solution 
for the 2BSDE with generator G (and thus also for the 2BSDE (|2.5p ) by passing to the 
limit in some sense in the sequence (Y n , Z n ). 

If we then follow the usual approach for standard BSDEs, the first step is to argue that 
thanks to the comparison theorem (which still holds true for Lipschitz 2BSDEs, see |24j). 
the sequence Y n is decreasing, and thanks to a priori estimates that it must converge 
Vh — Q-S- to some process Y. And this is exactly now that the situation becomes much more 
complicated with 2BSDEs. Indeed, if we were in the classical framework, this convergence 
of Y n together with the a priori estimates would be sufficient to prove the convergence in 
the usual H 2 space, thanks to the dominated convergence theorem. However, in our case, 
since the norms involve the supremum over a family of probability measures, this theorem 
can fail (we refer the reader to Section 2.6 in [23] for more details). Therefore, we cannot 
obtain directly that 



which is a crucial step in the approximation proof. 

This is precisely the major difficulty when considering the 2BSDE framework. The only 
monotone convergence Theorem in a similar setting has been proved by Denis, Hu and 
Peng (see [H]). However, one need to consider random variables X n which are regular in 
u), more precisely quasi-continuous, that is to say that for every e > 0, there exists an open 
set O e such that the X n are continuous in ui outside O s and such that 



Moreover, the set of probability measures considered must be weakly compact. This induces 
several fundamental problems when one tries to apply directly this Theorem to (Y n ) n >o- 

(i) First, if we assume that the terminal condition £ is in UC^ip), since the generator 
F (and thus G n ) are uniformly continuous in oj, we can reasonably expect to be able to 
prove that the Y n will be also continuous in uj, P — a.s., for every P € Vh- However, this 
is clearly not sufficient to obtain the quasi- continuity. Indeed, for each P, we would have 
a ¥ -negligible set outside of which the Y n are continuous in uj. But since the probability 
measures are mutually singular, this does not imply the existence of the open set of the 
definition of quasi- continuity. 

We moreover emphasize that it is a priori a very difficult problem to show the quasi- 
continuity of the solution of a 2BSDE, because by definition, it is defined P — a.s. for 
every P, and the quasi- continuity is by essence a notion related to the theory of capacities, 
not of probability measures. 

(ii) Next, it has been shown that if we assume that the matrices o p and if appearing in 
Definition \2.1\ are uniform in P, then the set Vh is on ^D weakly relatively compact. Then, 
we are left with two options. First, we can restrict ourselves to a closed subset of Vh, 
which will therefore be weakly compact. However, as pointed out in I25f . it is not possible 




sup P(O e ) < e. 
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to restrict arbitrarily the probability measures considered. Indeed, since the whole approach 
of \2J$ to prove existence of Lipschitz 2BSDEs relies on stochastic control and the dynamic 
programming equation, we need the set of processes a in the definition ofVs (that is to say 
our set of control processes) to be stable by concatenation and bifurcation (see for instance 
Remark 3.1 in J5fj in order to recover the results of [2^]. And it is not clear at all to us 
whether it is possible to find a closed subset of Vh satisfying this stability properties. 

Otherwise, we could work with the weak closure ofVn- The problem now is that the prob- 
ability measures in that closure no longer satisfy necessarily the martingale representation 
property and the 0-1 Blumenthal law. In that case (since the filtration ¥ will only be quasi- 
left continuous), and as already shown by El Karoui and Huang fl3f . we would need to 
redefine a solution of a 2BSDE by adding a martingale orthogonal to the canonical process. 
However, defining such solutions is a complicated problem outside of the scope of this paper. 

We hope to have convinced the reader that because of all the reasons listed above, it 
seems difficult in general to prove existence of a solution to a 2BSDE using approximation 
arguments. However, the situation is not hopeless. Indeed, in [23], the author uses such an 
approach to prove existence of a solution to a 2BSDE with a generator with linear growth 
satisfying some monotonicity condition. The idea is that in this case it is possible to show 
that the sequence of approximated generators converges uniformly in (y, z), and this allows 
to have a control on the difference \Y t n — Yt \ by a quantity which is regular enough to apply 
the monotone convergene Theorem of jllj . Nonetheless, this relies heavily on the type of 
approximation used and cannot a priori be extended to more general cases. 

Notwithstanding this, we will show an existence result in the next subsection using an 
approximation which has the particularity of being stationary, which immediately solves 
the convergence problems that we mentioned above. This approach is based on very recent 
results of Briand and Elie [8] on standard quadratic BSDEs. 

4.3 A stationary approximation 

For technical reasons that we will explain below, we will work throughout this subsection 
under a subset of Vh, which was first introduced in |26j . Namely, we will denote by S the 
set of processes a satisfying 

+oo +oo 

«*(<•>) = ^^a"' l l^(o;)l [TnHiTn+lH) (t), 

n=0 i=l 

where for each i and for each n, a n ' 1 is a bounded deterministic mapping, r n is an F-stopping 
time with To = 0, such that r n < r n+ i on {r n < +oo}, inf{n > 0, r n = +00} < +00, r n 
takes countably many values in some fixed Iq C [0, T] which is countable and dense in [0, T] 
and for each n, (Ef)i>\ C T Tn forms a partition of O. 

We will then consider the set Vh '■= £ Vh, a £ H) . As shown in [25], this set satisfies 
the right stability properties (already mentioned in the previous subsection) so much so 
that the Lipschitz theory of 2BSDEs still holds when we are working Vh — Q-s. Notice that 
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for the sake of simplicity, we will keep the same notations for the spaces considered under 
Vh or Vh- Let us now describe the Assumptions under which we will be working 

Assumption 4.1. Let Assumption \2.2\ holds, with the addition that the process (f> in (v) is 
bounded and that the mapping F is deterministic. 

The main result of this Section is then 

Theorem 4.1. Let Assumption \4-l\ hold. Assume further that £ € C^, that it is Malliavin 



differentiable Vh — Q-S- and that its Malliavin derivative is in B^. Then the 2BSDE (|2.5p 
(considered Vh — Q-s.) has a unique solution (Y, Z) 6 x M H . Moreover, the family 
{K ¥ , P £ Vh} can be aggregated. 

Proof. Uniqueness follows from Theorem 13. 1| so we concentrate on the existence part. Let 
us define the following sequence of generators 



F t n (y,z,a) :=F t [y 



\z\ A n 



z,a) , and F?{y,z) := F?(y,z,a t ) 



Then for each n, F n is uniformly Lipschitz in (y, z) and thanks to Assumption 14.11 we can 
apply the result of [23] to obtain the existence of a solution (Y n , Z n ) to the 2BSDE 



Y t n = i+ / F?(Y?,Z?)ds 



T 



Z"dB s + K^' n - K 



t ■> 



a.s., for all P G V H - (4.4) 



Moreover, we have for all P € Vh and for all t € [0, T] 
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P P.n 

ess sup y t , 

y ePi/(i+,P) 



a.s., 



(4.5) 



where (y ] 



\n JP,n\ 



is the unique solution of the Lipschitz BSDE with generator F n and 



terminal condition £ under P. Now, using Lemma 2.1 in [8] and its proof (see Remark 14.11 
below) under each P S Vh, we know that the sequence y p,n is actually stationary. Therefore, 
by (|4.5p . this also implies that the sequence Y n is stationary. Hence, we immediately have 
that Y n converges to some Y in DjJ. Moreover, we still have the representation 



Y, 



ess sup y t 

y 6^ H (t+,P) 



a.s., 



(4.6) 



Now, identifying the martingale parts in (|4.4|) . we also obtain that the sequence Z n is 
stationary and thus converges trivially in to some Z. For n large enough, we thus have 



Besides, we have by Assumption 14.11 



F t n (Yt,Z t ). 



F t n (Y t ,Z t 



< a + P\Y t \ + 



1 



l/2 \Zt\ An 



Z t 



Z t <a + /3\Y t \ + 



i 1/2 Z, 



, V H - q-s. 
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Since (Y, Z) S x M H , we can apply the dominated convergence theorem for the Lebesgue 
measure to obtain by continuity of F that 



F?(Y s n ,Z2)ds 



F s (Y s ,Z s )ds,V H ~ q-s. 



Using this result in (|4.4|) . this implies necessarily that for each P, K v,n converges P — a.s. 
to a non-decreasing process K p . Now, in order to verify that we indeed have obtained the 
solution, we need to check if the processes satisfy the minimum condition (|2.7p . Let 
P G Vh, t € [0, T] and P € V H (t A 
same notations 

-T 



) . From the proof of Theorem 13. 1| we have with the 



SYt. = Er 



M t dK[ 



ut 



t M T (M s )(K% -K\ 



Ef' 


£ ( Jo i^s + Vs)a s 1/2 dB s ^ 






£ 


I 




s )a s 1/2 dB s ^ 







For notational convenience, denote £t : = £ (^Jq ((j) s + r] s )a s l ^dB s j. Let r be the number 
given by Lemma 12.21 applied to £. Then we estimate 



<E? 



K 



^ inf (M s )(4 -Kf') 



l 

2r— 1 



<(SY t )^ < 



£ T \ 



r-l 
2r-l 



E 



E 



/ inf (M s )" 

t<s<± 



K 



inf (M s )~ 



2 

" r-l 



E 



K 



2r-l 



r-l 
2(2r-l) 



< C E: 



A" 



T 



r-l 
2(2r-l) 



(SYt) 



i 

2r-l 



By following the arguments of the proof of Theorem 13.11 (ii) and (iii) , we then deduce the 
minimum condition. Finally, the fact that the processes K F can be aggegated is a direct 
consequence of the general aggregation result of Theorem 5.1 in [26]. □ 



Remark 4.1. We emphasize that the result of Lemma 2.1 can only be applied when 

the generator is deterministic. However, even though F is indeed deterministic, F is not, 
because a is random. Nonetheless, given the particular form for the density of the quadratic 
variation of the canonical process we assumed in the definition of Vh, w & can apply the 
result of Briand and Elie between the stopping times and on each set of the partition ofQ, 
since then a and thus F is indeed deterministic. 



5 A pathwise proof of existence 

We have seen in the previous Section that it is usually extremely difficult to prove existence 
of a solution to a 2BSDE using monotone approximation techniques. Nonetheless, we have 
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shown in Theorem 13. II that if a solution exists, it will necessarily verify the representation 
(|2.7p . This gives us a natural candidate for the solution as a supremum of solutions to 
standard BSDEs. However, since those BSDEs are all defined on the support of mutually 
singular probability measures, it seems difficult to define such a supremum, because of the 
problems raised by the negligible sets. In order to overcome this, Soner, Touzi and Zhang 
proposed in [51] a pathwise construction of the solution to a 2BSDE. Let us describe briefly 
their strategy. 

The first step is to define pathwise the solution to a standard BSDE. For simplicity, let us 
consider first a BSDE with a generator equal to 0. Then, we know that the solution is given 
by the conditional expectation of the terminal condition. In order to define this solution 
pathwise, we can use the so-called regular conditional probability distribution (r.p.c.d. for 
short) of Stroock and Varadhan [27]. In the general case, the idea is similar and consists 
on defining BSDEs on a shifted canonical space. 

Finally, we have to prove measurability and regularity of the candidate solution thus ob- 
tained, and the decomposition (|2.5p is obtained through a non-linear Doob-Meyer decom- 
position. Our aim in this section is to extend this approach to the quadratic case. 

5.1 Notations 

For the convenience of the reader, we recall below some of the notations introduced in |24j . 

• For < t < T, denote by O* := {lo G C ([t,T],M d ) , w(t) = 0} the shifted canonical 
space, .B* the shifted canonical process, P the shifted Wiener measure and F* the filtration 
generated by B t . 

• For < s < t < T and uj G £l s , define the shifted path (J G ft 1 

w* := Lo r — oj-t, Vr G [t, T]. 

• For < s < t < T and u 6 !] s , w 6 Jl 1 define the concatenation path ui <S>t w G £l s by 

(w ® t u)(r) := w r l[ S)t) (r) + (uj t + a5 r )l[ ijT ](r), Vr <E [s,T\. 

• For < s < t < T and a J^-measurable random variable ^ on Q s , for each u> G £l s , define 
the shifted J^-measurable random variable ^* ,tJ on Or by 

f *•"(£;) :=£(uj® t uj), Vw G 0*. 

Similarly, for an F s -progressively measurable process X on [s,T] and (t,uj) G [s,T] x Q s , 
the shifted process |x* ,a; ,r G [t, T]| is F*-progressively measurable. 

• For a F-stopping time r, the r.c.p.d. of P (denoted P£) is a probability measure on Tt 
such that 

ffi?[f](w) =E p "[e], for P-a.e. w. 

It also induces naturally a probability measure P 7 "^ (that we also call the r.c.p.d. of P) 
on which in particular satisfies that for every bounded and J-r-measurable random 
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variable £ 

• We define similarly as in Section [5] the set "P^, by restricting to the shifted canonical 
space fi*, and its subset "P^. 

• Finally, we define our "shifted" generator 

F^{u,y,z) :=F s (w®i5,y,z,S*(5)), V(s,u5) G [t,T] xO*. 

Notice that thanks to Lemma 4.1 in [25], this generator coincides for P-a.e. uj with the 
shifted generator as defined above, that is to say 

F s (uo ® t uj,y,z,a s (uj ®tt5)). 

The advantage of the chosen "shifted" generator is that it inherits the uniform continuity 
in uj under the L°° norm of F. 

5.2 Existence when £ is in UCb(^) 

As mentioned at the beginning of the Section, we will need to prove some measurability and 
regularity on our candidate solution. For this purpose, we need to assume more regularity 
on the terminal condition. When £ is in UCb(^), by definition there exists a modulus of 
continuity function p for £ and F in uj. Then, for any < t < s < T, (y, z) G [0, T\ x M x R rf 
and uj,uj' G J), uj £ f2*, 

< p(||u; -u/y and |f*>" (w.y.z) - F^' (uj,y,z)\ < p (\\uj - uj'\\ t ) , 
where ||w|| t := sup 0<s<t \uj s \ , < t < T. 

To prove existence, as in [23], we define the following value process Vt pathwise 

V t (uj) := sup yf AuJ (T,0 , for all (t,uj) E [0,T] x n, (5.1) 

where, for any (t^uj) G [0,T] xQ, F G i 2 G [ti,T], and any ^-measurable 7] G L°° (P), 
we denote y t ^ tl,u {t2,rj) := ytf 1,0J , where (y p,tl ' u , ;z p >* 1 > a ') is the solution of the following 
BSDE on the shifted space il* 1 under P 



/* ^2 

^>*i><^ - y z y^dBl\ se[h,t 2 },¥- a.s. (5.2) 



We recall that since the Blumenthal zero-one law holds for all our probability measures, 
yf' t,Ld (1,£) is constant for any given (t, uj) and P G Vjj. Therefore, the process V is well 
defined. However, we still do not know anything about its measurability. The following 
Lemma answers this question and explains the uniform continuity Assumptions in uj we 
made. 
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Lemma 5.1. Let £ 6 UCb(f2). Under Assumption \2.1\ or Assumption [2TM with the addition 
that the -norms of £ and F° are small enough, we have 

|T4H|<c(l+||e|| LS ), for all (t, oj) G [0, T] x ft. 

Furthermore, \V t (a/) - Vj (u/)| < Cp(\\u - u)'\\ t ) , for all (t,uj,u)') G [0,T] x fl 2 . In partic- 
ular, V t is T t -measurable for every t G [0,T]. 



Proof, (i) For each (t, oj) G [0, T]xfl and P G 7^, note that for s G [i, T], we have P - a.s. 



e 



T r 



(0) + XrVr + (Vr + <M a 



\ 1/2 P,t, ( 



where A is bounded and r] satisfies 

\Vr\ < V 



-1/2 P,i,cj 



a.s. 



Then proceeding exactly as in the second step of the proof of Theorem 13.11 we can define 
a bounded process M and a probability measure Q equivalent to P such that 



lit 



<E^ [m t \^\) <c(i + u\ 



By arbitrariness of P, we get |Vt(w)| < C(l + ||£||l<x>). 

(ii) The proof is exactly the same as above, except that we need to use uniform continuity 
in u of £ t ' w and F*' w . In fact, if we define for {t,uj,uj') G [0,T] x fl 2 



',t,u 



5y := - tfW , 5z := z l - Z' 



then we get with the same notations 



''*> w ' ■= — ^*' a;/ SF ■= F t,u) — F t,UJ ' 



\8yt\=W> 



M T 5i+ / M s 8F s ds 



<Cp(\\u-oj'\\ t ). 



We get the result by arbitrariness of P. 

Then, we show the same dynamic programming principle as Proposition 4.7 in [25 



□ 



Proposition 5.1. Let £ G UCb(O). Under Assumption \2.1\ or Assumption \2.2\ with the 
addition that the -norms of £ and F° are small enough, we have for all < t\ < ti < T 
and for all oj G ft 

The proof is almost the same as the proof in [25] , but we give it for the convenience of the 
reader. 

Proof. Without loss of generality, we can assume that t\ = and ti = t. Thus, we have 
to prove 

v (u) = sup y$(t,v t ). 
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Denote := Q?(T, £), Z P (T, £)) 



(i) For any P G Pff, it follows from Lemma 4.3 in |25| . that for P — a.e. uj G fi, the r.c.p.d. 
P*> w eP^. By Tevzadze [28], we know that when the terminal condition and F are small, 
quadratic BSDEs whose generator satisfies Assumption (|2.2p can be constructed via Picard 
iteration. Thus, it means that at each step of the iteration, the solution can be formulated 
as a conditional expectation under P. Then, for general terminal conditions, Tevzadze 
showed that if the generator satisfies Assumption ()2. 1|) (v), the solution of the quadratic 
BSDE can be written as a sum of quadratic BSDEs with small terminal condtion. By the 
properties of the r.p.c.d., this implies that 

y f (u) = yf (T,0, for P - a.e. uj G 0. 

By definition of Vt and the comparison principle for quadratic BSDEs, we deduce that 
Uq < 3^o Vt) an d it follows from the arbitrariness of P that 

V (u) < sup 3^o (t,V t ). 

f&V H 

(ii) For the other inequality, we proceed as in [25]. Let P G Vh an d e > 0. The idea is to 
use the definition of V as a supremum to obtain an e-optimizer. However, since V depends 
obviously on uj, we have to find a way to control its dependence in uj by restricting it in a 
small ball. But, since the canonical space is separable, this is easy. Indeed, there exists a 
partition {E£)i>\ C Ft such that \\uj — uj'\\ t < e for any i and any uj,uj' G E\. 

Now for each i, fix an G El and let, as advocated above, P£ be an e— optimizer of Vt(uji). 
If we define for each n > 1, P n := P™' 6 by 



F n (E) :=E ff 



i=l 



t.co 



+ F(E n E?), where E? := \J l>n E\ 



t • 



then, by the proof of Proposition 4.7 in [25], we know that P n G Vh and that 

V t < yT + e + Cp(e), P n - a.s. on U™ =1 E\. 



Let now (y n ,z n ) := (y n,e ,z n,e ) be the solution of the following BSDE on [0, t] 



yf+e + Cpie) 



Z^dB r 



a.s. 



(5.3) 



Note that since P n = P on Ft, the equality (|5.3p also holds P — a.s. By the comparison 
theorem, we know that y$ (t, Vt) < yfi . Using the same arguments and notations as in the 
proof of Lemma l5.lt we obtain 



e + p{?) + 



Vt 



Then, by Lemma 1 5. H we have 

3^ V t ) < y'S < V (uj) + C (e + p(e) + E ( 
The result follows from letting n go to +oo and e to 0. 



1, 



M, 



□ 
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Remark 5.1. We want to emphasize here that it is only because of this Proposition prov- 
ing the dynamic programming equation that we had to consider Tevzadze J28j approach to 
quadratic BSDEs, instead of the more classical approach of Kobylanski [PT\[ . Indeed, as 
pointed out in the proof, for technical reasons we want to be able to construct solutions of 
BSDEs via Picard iterations, to build upon the known properties of the r.c.p.d. Using the 
Assumptions ]^. 1\ or \2.'2\ with the addition that the -norms of £ and F° are small enough, 
this allows us to recover this property. 

Now that we solved the measurability issues for Vt, we need to study its regularity in time. 
However, it seems difficult to obtain a result directly, given the definition of V. This is the 
reason why we define now for all (t,u), the F + -progressively measurable process 

V t + := Ihn V r . 

reQn(t,T],rit 

This new value process will then be proved to be cadlag. Notice that a priori V + is only 
F + -progressively measurable, and not F-progressively measurable. This explains why in the 
definition of the spaces in Section 12.3. 1\ the processes are assumed to be F + -progressively 
measurable. 

Lemma 5.2. Under the conditions of the previous Proposition, we have 

V + = lim V r , Vn-q.s. 

reQn(t,T],rlt 

and thus V + is cadlag Vh — qs. 

Proof. Actually, we can proceed exactly as in the proof of Lemma 4.8 in [25], since the 
theory of ^-expectations of Peng has been extended by Ma and Yao in [18] to the quadratic 
case (see in particular their Corollary 5.6 for our purpose). □ 



Finally, proceeding exactly as in Steps 1 and 2 of the proof of Theorem 4.5 in [25] . and in 
particular using the Doob-Meyer decomposition proved in [18J (Theorem 5.2), we can get the 
existence of a universal process Z and a family of non-decreasing processes {K v ,¥ £ Vh} 
such that 

v t + = V + + / F*(V+, Z s )ds + / Z s dB s - Kf, P - a.s. VP € V H - 
Jo Jo 

For the sake of completeness, we provide the representation (|3.2p for V and V + , and that, 
as shown in Proposition 4.11 of [25] . we actually have V = V + , Vh — Q-s., which shows 
that in the case of a terminal condition in UC{,(£1), the solution of the 2BSDE is actually 
F-progressively measurable. This will be important in Section [7] 

Proposition 5.2. Let £ € UCb(£l). Under Assumption \2.1\ or Assumption \2.2\ with the 
addition that the -norms of ^ and F° are small enough, we have 

V t = esssup p yf'{T,£) and V t + = esssup p $ (T,£), ¥ - a.s., VP G Vh- 

v'eVH(t,¥) p'e-Pff(t+,P) 

Besides, we also have for all t, Vt = V t + , Vh — Q-s. 
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Proof. The proof for the representations is the same as the proof of proposition 4.10 in 
[25], since we also have a stability result for quadratic BSDEs under our assumptions. For 
the equality between V and V + , we also refer to the proof of Proposition 4.11 in |25j . □ 

To be sure that we have found a solution to our 2BSDE, it remains to check that the family 
of non-decreasing processes above satisfies the minimum condition. Let P G Vh, t G [0, T] 
and P G Vh (t + , P) • From the proof of Theorem 13.11 we have with the same notations 



SV t = E 



MtdKf 



t mf T (M s )(^ 



K 









)a~ 1/2 dB s ^ 


t< inyM s )(^'-Af ) 


Ef 









s )a s l/2 dB s ^ 







We can proceed exactly as in the proof of Theorem 14.11 in Subsection 14.3 



□ 



Remark 5.2. In order to prove the minimum condition it is fundamental that the process 
M above is bounded from below. For instance, it would not be the case if we had replaced 
the Lipschitz assumption on y by a monotonicity condition as in \2c 



5.3 Main result 

We are now in position to state the main result of this section 

Theorem 5.1. Let £ G C^. Under Assumption \2. 1\ or Assumption \2.2\ with the addition 
that the -norms of £ and F° are small enough, there exists a unique solution (Y,Z) G 
x W 2 H of the 2BSDE rf23j) . 

Proof. For £ G £<S, there exists £ n G UC b (ft) such that ||f - £ n || -> 0. Then, thanks 

n— >+oo 

to the a priori estimates obtained in Proposition [372], we can proceed exactly as in the proof 
of Theorem 4.6 (ii) in [24] to obtain the solution as a limit of the solution of the 2BSDE 
(I2.5P with terminal condition £ n . □ 



6 An application to robust risk-sensitive control 

One application of classical quadratic BSDEs is to study risk-sensitive control problems, 
see El Karoui, Hamadene et Matoussi |14] for more details. In this section, we will consider 
a robust version of these problems. 

First of all, for technical reasons, we restrict the probability measures in Vh '■= VsCWh, 
where Vs is defined in Subsection 12.11 Then a is uniformly bounded by some a, a G 

For each P G Vh, we can define a P-Brownian motion W v by 

dWf = a; 1/2 dB t P - a.s. 

Let us now consider some system, whose evolution is decribed (for simplicity) by the canon- 
ical process B. A controller then intervenes on the system via an adapted stochastic process 
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u which takes its values in a compact metric space U. The set of those controls is called ad- 
missible and denoted by U. When the controller acts with u under the probability P E Vh, 
the dynamic of the controlled system remains the same, but now under the probability 
measure P" defined by its density with respect to P 



~d¥ 



T 



-1/2 - 



exp [ I a t g(t,B.,u t )dW; - - 
lo z Jo 



T -1/2 

a t g(t,B.,u t 



dt 



where g(t,u),u) is assumed to be bounded, continuous with respect to u, adapted and 
uniformly continuous in u. Notice that this probability measure is well defined since a is 
uniformly bounded. 

Then, under ¥ u , the dynamic of the system is given by 



dB t = g(t, B., u t )dt + a t ' dW 
where W^ ,u is a Brownian motion under P" defined by 



a.s. 



dwf' u = dWf 



-1/2 



g(t,B.,u t )dt. 



When the controller is risk seeking, we assume that the reward functional of the control 
action is given by the following expression 



Vn G U, J{u) := sup E 



exp^ej h(s,B.,u s )ds + ^(B T ) 



where 6 > is a real parameter which represents the sensitiveness of the controller with 
respect to risk. Here h(t,u),u) is assumed to adapted and continuous in u, and both ^ 
and h are assumed to be bounded and uniformly continuous in u for the ||-|| norm. We 
are interested in finding an admissible control u* which maximizes the reward J(u) for the 
controller. 

We begin with establishing the link between J(u) and 2BSDEs in the following proposition 

Proposition 6.1. There exists a unique solution (Y U ,Z U ) of the 2BSDE associated with 
the generator zg(t, B.,ut) + h(t, B.,ut) + ||a^ 2 z| 2 , i.e., P — a.s., for all P G Vh 

Y t u = V(B T ) + ^ (z»g(s,B.,u s ) + h(s,B.,u s ) + ^a^Z^ ds - £ Z u s dB s - dKf '. 

(6.1) 

Moreover J(u) = exp (6Yq). 



Proof. With our assumptions on g, h and we know that the generator satisfies the 
Assumption 12 . 1 [ therefore there exists a unique solution to the 2BSDE (|6.ip . According to 
|14j . the solution to the classical BSDE with the same terminal condition and generator as 
the 2BSDE JO) under each P is 



Vt = z In I Ei 



exp^ej^ h(s,B.,u s )ds + W(B T ) 



a.s. 



2S 



Then by the representation for Y u , we have 



expfej^ h(s,B.,u s )ds + *f>(B T ) 



Yt= -, esssup r In E^ 

Since the functional ln(x) is monotone non-decreasing, then 
! / r / rT 



a.s. 



Y t u = - In esssup r Ef 



a.s. 



expfe f h(s,B.,u s )ds + ^{B T ) 
v -'ePH(*+,P) L ^ ^ 

Therefore, we have J(it) = exp {^Yq 4 }. □ 

As explained in p3], by applying Benes' selection theorem, there exists a measurable version 
u*(t,B.,z) of 

arg max I(t, B.,z, u) := zg(t, B.,u) + h(t, B., u). 

We know that I*(t,B.,z) := sup u£U I(t, B., z, u) = I(t, B., z,u*(t, B., z)) is convex uni- 
formly Lipschitz in z because it is the supremum of functions which are linear in z. So the 
mapping z — > I*(t, B., z) + ^\cL t ^ 2 z\ 2 is continuous with quadratic growth, implying that a 
solution (y*' F ', z*^) of the BSDE associated to this generator exists. Then we have 

Theorem 6.1. There exists a unique solution (Y*,Z*) to the following 2BSDE 

Y* = y(B T ) + jT (l*(s, B., Z* s ) + ^|ay 2 Z s *| 2 ^ ds - £ Z*dB s + K*/ - Kf. (6.2) 

The admissible control u* := (u*(t,B.,Z^))t<T is optimal and (exp(Y t *))t<T is the value 
function of the robust risk-sensitive control problem, i.e., for any t <T we have: 



exp(Y 4 j = ess sup ess sup E t 

P'ePjj(t+,P) u&a 



exp [9 j h(s,B.,u s )ds + $>(B T ) 



Proof. First, we need to prove the existence of a solution to the quadratic 2BSDE (|6.2p . 
Unlike in Proposition 16. 1\ here u* also depends on z, so we do not know whether I* is twice 
differentiable with respect to z. Therefore the generator of the 2BSDE may not satisfy the 
Assumption 12.11 But it's easy to see that it always satisfies the weaker Assumption 12.21 
and we only need this Assumption to have uniqueness of the solution. Moreover, it was 
also the only one used to prove the minimum condition for the familly of non-decreasing 
processes in Subsection 15.21 Therefore, exactly as in Section [H for P G Vh, by making the 
exponential change 



Yt 



BY* -79 

e * , Z t 



9Y t Z* t , K t ■, 



9 / Y s dKl 



o 



0<s<t 



JY* 



QY* 

e b- 



9AY*e eY s* 



we see that (Y, Z,K ) formally verifies the following equation 

pT pT 

y i = e e*(B T ) + / sup {z s g{s, B., u) + 9Y s h(s, B., u)} ds- / ~Z s dB s +Tf T -lf t , F-a.s. 

(6.3) 
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Since this is 2BSDE with Lipschitz generator from Soner, Touzi and Zhang [24j . we know 
that (Y,Z, K ) exists, is unique and satisfies the representation property (|3.2p . Arguing 
exactly as in Subsection 14.11 for the purely quadratic 2BSDEs, we can then obtain the 
existence. Now, from p3], we have that 

exp (y*' P ) = esssup p Ef exp (o h(s, B., u s )ds + ^(B T )^j . 

Then the representation for Y* implies the desired result. □ 



7 Connection with fully nonlinear PDEs 

In this section, we assume that all the randomness in H only depends on the current value 
of the canonical process B (the so-called Markov property) 

H t {w,y,z,i) = h(t,B t (u),y,z,-y), 

where h : [0, T] x M. d x R x R rf x — > R is a deterministic map. Then, we define as in 
Section [5] the corresponding conjugate and bi-conjugate functions 

f(t, x, y, z, a) := sup { -Tr [07] - h(t, x, y, z, 7) 1 (7.1) 
■yeD h 12 J 

h{t,x,y,z,j) := sup { -Tr [aj] - f(t, x, y, z, a) \ (7.2) 
aes>° I 2 J 

We denote Vh '■= Vh, and following [23] . we strengthen Assumptions 12.11 and 12.21 

Assumption 7.1. (i) The domain Df t of f in a is independent of (x,y,z). 

(ii) On Df t , f is uniformly continuous in t, uniformly in a. 

(iii) / is continuous in z and there exists (a, f3, 7) such that 

\f(t,x,y,z,a)\ < a + P\y\ + |ja 1/2 z| 2 , for all (t,x,y,z,a). 

(iv) / is uniformly continuous in x, uniformly in (t,y,z,a), with a modulus of continuity 
p which has polynomial growth. 

(v) / is C 1 in y and C 2 in z, and there are constants r and 9 such that for all (t, x, y, z, a) 

\D y f(t, x, y, z,a)\ < r, \D z f(t, x,y,z,a)\ < r + 6\a 1/2 z\, \D 2 zz f(t, x, y, z,a)\ < 9. 
Assumption 7.2. Let points (i) through (iv) hold and 

(v) There exists [i > and a bounded M. d -valued function (p such that for all (t, y, z, z' , a) 

\f(t,x,y,z,a)-f(t,x : y : z',a)-0(t).a 1/2 (z-z')\ < fia 1/2 \z - z'\ ( a l/2 z + a l/2 z' ) . 

(vi) / is Lipschitz in y, uniformly in (t,x,z,a). 
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Let now g : R d -> R be a Lebesgue measurable and bounded function. Our object of 
interest here is the following 2BSDE with terminal condition £ = g(Bx) 

Y t =g{B T )+ [ T f(s,B s ,Y s ,Z s ,al /2 )ds- C Z s dB s + - K? ', V h -q.s. (7.3) 
Jt Jt 



The aim of this section is to generalize the results of [21] and obtain the connection Yt = 
v(t,Bt), Vh — Q-S-, where v verifies in some sense the following fully nonlinear PDE 

§2(t, x) + h (t, x, v{t, x),Dv{t, x),D 2 v(t, x)) = 0, t G [0, T) 

(7.4) 

v(T,x) = g(x). 

Following the classical terminology in the BSDE litterature, we say that the solution of 
the 2BSDE is Markovian if it can be represented by a deterministic function of t and B t . 
In this subsection, we will construct such a function following the same spirit as in the 
construction in the previous section. We want to emphasize that the proofs here follow 
very closely the proofs in [21], so we will sometimes only sketch them. 

With the same notations for shifted spaces, we define for any (t, x) G [0, T] x W d 

BY :=x + B l s , for all s G [t,T]. 

Let now r be an F*-stopping time, P G V\ and rj a P-bounded J^-measurable random 
variable. Similarly as in (j5.2|) . we denote (y p ' t ' x , z p,t > x ) := (y p ' t ' x (T,rj),Z p < t ' x (r,r])) the 
unique solution of the following BSDE 



+ [ ' f{u,BY,^ x X' t>x ,K)du- f zl^dB^, t<s<r, F-a.s. (7.5) 

J s J s 

Next, we define the following deterministic function (by virtue of the Blumenthal — 1 law) 
u(t,x) := sup3 ; t P '*' X (r,5(^ x )), for (t,x) G [0,T] x R d . (7.6) 

We then have the following Theorem, which is actually Theorem 5.9 of |24j in our framework 



Theorem 7.1. Let Assumption 7.1 hold, and assume that g is bounded and uniformly 
continuous. Then the 2BSDE (|7.3p has a unique solution (Y, Z) G x and we have 
Yt = u(t, Bt). Moreover, u is uniformly continuous in x, uniformly in t and right- continuous 
in t. 

Proof. The existence and uniqueness for the 2BSDE follows directly from Theorem 15.11 
Since £ G UCb(^), we have with the notations of the previous section V t = u{t,B t ). But, 
by Proposition 15.21 we know that Yt = Vt, hence the first result. 

Then the uniform continuity of u is a simple consequence of Lemma 15. 1[ Finally, the 
right-continuity of u in t can be obtained exactly as in the proof of Theorem 5.9 in [24J. □ 
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7.1 Non-linear Feynman-Kac formula in the quadratic case 

Exactly as in the classical case and as in Theorem 5.3 in [24] , we have a non-linear version 
of the Feynaman-Kac formula. The proof is the same as in [25], so we omit it. Notice 
however that it is more involved than in the classical case, mainly due to the technicalities 
introduced by the quasi-sure framework. 



Theorem 7.2. Under Assumption \7.1\ suppose that h is continuous in its domain, that Df 
is independent oft and is bounded both from above and away from 0. Let v € C 1,2 ([0, T), ] 
be a classical solution of (|7.4|) with {(v, Dv)(t, B t )} 0<t<T £ x H^. Then 



Y t := v(t,B t ), Z t := Dv(t,B t ), K t := I 

Jo 



Jxicds . 



is the unique solution of the quadratic 2BSDE (|7.3p . where 



h := h(t, B U Y U Z u T t ) - ^Tr [a t 1/2 r t ] + f(t, B t ,Y t , Z t ,a t ) and T t := D 2 v(t, B t ). 
7.2 The viscosity solution property 

As usual when dealing with possibly discontinuous viscosity solutions, we introduce the 
following upper and lower-semicontinuous envelopes 

u*(t,x) := lim u(t',x'), u*(t,x) := lim u(t',x') 

(t' ,x')->(t,x) (f ,x')->{t,x) 

%*(•&):= lim h(tf),h*(#):= li h(#) 

In order to prove the main Theorem of this subsection, we will need the following Proposi- 
tion, whose proof (which is rather technical) is omitted, since it is exactly the same as the 
proof of Propositions 5.10 and 5.14 and Lemma 6.2 in 



Proposition 7.1. Let Assumption ^ .1\ hold. Then for any bounded function g and (t,x) 
(i) V {t p ,P € Vl), ^-stopping times , we have u(t,x) < sup yf '*' x (t p , u*(t f , B^)). 



(ii) If in addition g is lower-semicontinuous, then u(t,x) = sup y t ' ' x (t ,u(t ,B '£)). 

Now we can state the main Theorem of this section 
Theorem 7.3. Let Assumption \7. l\ hold true. Then 

(i) u is a viscosity subsolution of 

-d t u* -h*(-,u*,Du*,D 2 u*) < 0, on [0,T) x R d . 

(ii) If in addition g is lower-semicontinuous and Df is independent of t, then u is a 
viscosity supersolution of 

-dtu* -h t (;u„ Du„ D 2 u*) > 0, on [0, T) x R d . 
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Proof. The proof follows closely the proof of Theorem 5.11 in [23], with some minor 
modifications (notably when we prove (|7.10p ). We provide the proof of the supersolution 
property for the convenience of the reader, the subsolution property can be proved similarly. 
Assume to the contrary that for some (to,xo) € [0, T) x M d , 

= (u* - fifa, x ) > (u* - <P)(t, x) for all (t, x) € [0, T) x R d \ {(to, x )} , (7.7) 

(-d t <t>-h*{;<i>,D<l>,D 2 <t>j) (to,x o )>0, (7.8) 

for some smooth and bounded function 4> (we can assume w.l.o.g. that (f> is bounded since 
we are working with bounded solutions of 2BSDEs). Now since <fi is smooth and since by 
definition h* is upper-semicontinuous, there exists an open ball 0(r, (to, xq)) centered at 
(to,xo) with radius r, which can be chosen less than T — to, such that 

-8t<f> - h(-, 0, D<f>, D 2 cp) > 0, on 0(r, (t , x )). 

By definition of h, this implies that for any a € S^ -0 

- dt(f> - ^Tr [aD 2 <f>] + /(•, 0, D<f>, a) > 0, on 0(r, (t ,x )). (7.9) 
Let us now denote a := — max (u* — </>). By (17. 71) . this quantity is strictly positive. Let 

dO(r,(t ,x )) 

now (t n ,x n ) be a sequence in O(r,(to,xo)) such that (t n ,x n ) — > (t ,xo) and u(t n ,x n ) — > 
u*(to,xo). Denote the following stopping time 

r n :=inf{s>i n , (s,Bl n ' Xn £ 0(r, (t , x ))} . 

Since r < T — to, we have r n < T and thus (r n , B T "' Xn ) € dO(r, (to,xo)). Hence, we have 

c n := (<j)-u)(t n ,x n ) ^0 and u*(r n , B^' Xn ) < 0(r n , B%< Xn ) - fi. 

Fix now some P n € V l £ . By the comparison Theorem for quadratic BSDEs, we have 

yf;^(r n ,u*(T n ,Bi : ^)) < y^' tn ' Xn (r n , <f>(r n , B%f n ) - n). 

Then proceeding exactly as in the second step of the proof of Theorem 13.11 we can define 
a bounded process M n , whose bounds only depend on T and the Lipchitz constant of / in 
y, and a probability measure (Q^ equivalent to IP72 such that 

^ P n "' ta (r n ,0(r n)J B^»)- M ) -y^^rnA^B^)) = -E^[M Tn »] < -J, 

for some strictly positive constant fi' which is independent of n. Hence, we obtain by 
definition of c n 

yf:' tn ' Xn (r n ,u*(r n ,Bl^)) - u(t n ,x n ) < yf™' tn ' Xn (T n , <f>(r n , B^ Xn )) - <f>(t n ,x n ) + c n - . 

(7.10) 

With the same arguments as above, it is then easy to show with Ito's formula that 



y^ Xn (T n ,<t>(T n ,B^))-^(t n ,x n ) = ^\ 



" M^ n s ds 



33 



where C := (-d t 0-±Tr [a^D 2 ^] +f(-,D<j), a*))(s, But by £Ll and the definition 

of r n , we know that for t n < s < r n , tp™ > 0. Recalling fj7. lOf) . we then get 

^^(r n ,u m {T n ,B^)) - u(t n ,x n ) <c n - n'. 

Since c n does not depend on P n , and since the right-hand side is strictly negative for n 
large enough, this contradicts Proposition 17. If i). □ 
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